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Abstract. In |17| we obtained a nonconventional invariance principle (func- 
tional central limit theorem) for sufficiently fast mixing stochastic processes 
with discrete and continuous time. In this paper we derive a nonconventional 
invariance principle for sufficiently well mixing random fields. 



1. Introduction 

Nonconventional ergodic theorems (see |!l2]) known also after [3] as poli- 
nomial ergodic theorems studied the limits of expressions having the form 

where T is a weakly mixing measure preserving 
transformation, /^'s are bounded measurable functions and g^'s are polynomials 
taking on integer values on the integers. Originally, these results were motivated 
by applications to multiple recurrence for dynamical systems taking functions fi 
being indicators of some measurable sets. Later such results were extended to the 
case when q^'s are polinomials on Z'' (see 18 ) and to some actions (see |^). 

Using the language of probability this kind of results may be called nonconven- 
tional laws of large numbers and as a natural follow up we arrived at the invariance 
principle (functional central limit theorem) in |17| showing convergence in distri- 
bution to Gaussian processes for expressions of the form 

(1.1) 1/Vn {F{X{q,{n)),...,X{q,{n)))~F) 

0<n<Nt 

where X(n), n > is a sufficiently fast a, p or i/j-mixing vector valued process with 
some moment conditions and stationarity properties, F is a continuous function 
with polinomial growth and certain regularity properties, F = J Fd{^i x ■■ ■ x fi), 
fi is the distribution of each X (n) , qj (n) = jn, j < k < £ and qj, j = k + 1, ...,£ are 
positive functions taking on integer values on integers with some growth conditions 
which are satisfied, for instance, when g^'s are polynomials of growing degrees. 

The goal of this paper is to prove an invariance principle type result when n G Z"^ 
is multidimensional. This can be done either by considering functions : Z'' -> Z-)_ 
with X{n),n > being again a vector valued stochastic process or, more generally, 
considering maps qi : Z'' — > Z'^ taking now X{n),n G Z'^ to be a vector valued 
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random field which will be our setup in this paper. Namely, for t — (ti, ...,tu) e 
[0, 1]" and a positive integer N we consider expressions of the form 

(1.2) ^^(t) = 7V-^/2 {F{X{q^{n)),...,X{q,{n)))~F) 

n={ni,...,n^):0<ni<Nti Vi 

where X{n),n € is a sufficiently well mixing vector valued random field, with 
some moment conditions and stationarity properties, F and F are similar to above, 
qj{n) = jn, j < k and qi'.l/ ^ Z", j = fc + 1, map Z"^ — {n = (ni, ...,n^) € 
Z" : > 0, i = 1, i^} into itself. Assuming some growth conditions of % > k 
in \n\ we will show that the random field ^Ar(t) converges in distribution to a 
Gaussian random field on [0, 1]"^. 

In [17] we were able to obtain the latter result for one dimensional n relying on 
martingale approximations and martingale limit theorems but for random fields this 
machinery is not readily available. Still, we are able to combine some of mixingale 
technique from [19] and [20j together with an appropriate grouping of summands 
in (jl.2p in order to obtain both convergence of finite dimensional distributions 
and the tightness of infinite dimensional ones. Other known methods which work 
successfully when proving limit theorems for random fields (see, for instance, [S], [7], 
[8j and [21] ) rely one way or another on characteristic functions (or other devices 
based on weak dependence) which are hard to deal with in the nonconventional 
setup as demonstrated in |15j in view of the strong dependence of the summands 
in (|1.2p on the far away members of the random field. For specific lattice models 
with sufficiently good mixing properties to fit our setup we refer the reader to [T] 
and references there. 

2. Preliminaries and main results 

Our setup consists of a p-dimensional vector random field {X(n), n e 
, X{n) e M^} on a probability space {fl,J-,P) and of a family of cr-algebras 
J"r C J", r C Z" such that J"r C J'a if T C A C Z". It is often convenient to 
measure the dependence between two sub cr-algebras Q^T-L <Z F via the quantities 

(2.1) rUq^p{Q, %) ~ sup{||i?(g|t/) — Eg\\p : g is V. — measurable and ||,g||q < 1} 

where the suprcmum is taken over real functions and || • || r is the L^{^1, F, P)-norm. 
Then more familiar a, p, 4> and V'-niixing (dependence) coefficients can be expressed 
via the formulas (see [6J , Ch. 4 ) , 

(t>iG,H) = \tiJoo,oo{Q,'H) and ^{G.'H) = ^\,oo{Q,'H). 

We set also 

(2.2) w,,p(r)= sup (|ru Ariro,,p(J-r,J"A)) 

r,A: dist(r,A)>i- 

where F and A are finite nonempty subsets of Z"^, ciist(F, A) — inf„gr,7T,GA \n — n\ 
and we write |F| for cardinality of a set F while, as usual, for numbers or vectors 
I • I will denote their absolute values or lengths. As shown in |10j imposing decay 
conditions on dependence coefficients which do not take into account sizes of sets F 
and A as in ()2.2p would exclude from our setup simple examples of Gibbs random 
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fields. Define also 

"(0 = ^^oo4(0: P(0 = ^2,2(1), 0(0 = 1^00,00(0 and ^^(0 = Tni^oo(0- 
Our setup includes also conditions on the approximation rate 

(2.3) /3p(r) = sup \\X{n) - E{Xin)\Tu^^„^)\\p 

where Ur{n) = {h G Z'^ : \n — n\ < r} is the r-neghborhood on n in Z"^. Further- 
more, we do not require stationarity of the random field X{n),ji G Z'' assuming 
only that the distribution of X{n) does not depend on n and the joint distribution 
of {X{n), X(n')} depends only on n — n' which we write for further references by 

(2.4) X{n) - ^ and (X(n), X{n')) - /x„_„' for ah n, n 
where Y Z means that Y and Z have the same distribution. 



Next, let F = F(xi, cc^), Xj € R*^ be a function on MP such that for some 
i, if > 0, K e (0, 1] and all cc^, t/i e M*^, « = 1, ^, we have 

I e e 

(2.5) |F(xi, ...,x,) - F(yi, y,)| < if (l + ^ + ^ ly.f) ^ \xj - y,r 

i=i i=i j=i 

and 



(2.6) |F(xi,...,Xf)| <if(l- 



E 



Our assumptions on F enable us to include, for instance, products F{xi, ...,Xi) — 
2^112:22 • • • Xii, where Xi = {xn, xn) G K^, which is sometimes useful. To simplify 
formulas we assume a centering condition 

(2.7) F = J F{xi,...,xi)dn{xi)---dfi{xe)^0 

which is not really a restriction since we always can replace F hy F — F. 

Our goal is to prove an invariance principle (functional central limit theorem) 
for £,Nit), t e [0, 1]" defined by pr2)) where qj{n) = jn if j = 1, 2, k < £ for some 
given positive integers k, £ and if fc < £ then qj : Z'^ 'Z'^ , j = k + 1, ...,£ satisfy the 
conditions below. Before we formulate them observe that already the case k = £ is 
of major interest and we add q'jS with j > k which grow faster than linearly mostly 
for the sake of completeness which under appropriate assumptions does not cause 
substantial problems. We assume that |(7i(n)| < |g2(?^)| < ••■ < |'?£(?^)| whenever 
\n\ ^ and n G Z"^ = {m = (toi, m,y) G U' : nii > OVi}. Furthermore, we 
assume that for fc + 1 < i < 

(2.8) inf (|g,(n)|-|g,(n)|)(|n|-H)-i >0, 

\n\>\n\ 

(2.9) lim inf (\qi{h) ~ qi{n)\ — \n — n\) — 00, 

n|— >oo h:h^n 

(2.10) qi{n) ^ qi{h) if n ^ n, lim min [\qi{n)\ - \qi{n)\ - \n\) ^ 00 

|n|— J-oo l<i 

and for any e > 0, 

(2.11) lim inf min (|(7i(n)| — |(7i(n)| — |n — n|) = 00. 

|n|— )-oo n: |fi|>e|n| l<i 
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For each > set 

(2.12) ^l^\\X\\l = E\X[n)\' ^ J N%. 
Our main result relies on 

2.1. Assumption. With d — {£ — l)p there exist p,q > I, m > A and S > with 
S < K, pK > d satisfying 

oo 

(2.13) ^?5''z:j7,,p(0=%,g)<oo, 

1=0 

oo 

(2.14) ^r5-/3^'(r) <oo, 

. 1 1 i + 1 (5 

(2.15) 7m < oo, 72„i < oo with - > - H h -. 

2 p 771 g 

In order to give a detailed statement of our main result as well as for its proof 
it will be essential to represent the function F = F{xi, X2, ■ ■ ■ , xe) in the form 

(2.16) F = Fi{xi) + ■■■+ Fe{xi,X2, . . . ,x,) 
where for i < £, 

(2.17) Fi{xi,...,Xi) = J F{xi,X2, . ■ . ,xe) dfi{xi+i) ■ ■ ■ dfi{xi) 

- J F{xi,X2, ■ --^xe) d^{xi) ■ ■ ■ dfi{xi) 

and 

Fi{xi,X2, ■ ■ ■ ,xe) = F{xi,X2,...,Xi) - J F{xi,X2, ■ ■ ■ ,Xi) dij,{xi) 
which ensures, in particular, that 

(2.18) J Fi{xi,X2,...,Xi-i,Xi)d^{xi) = V .ti, X2, . . . , 

We write t — (ii, ...jt^) > s = (si, s^) if ti > Si for all i and for such s,t E [0, 1]'' 
we set A]\f{s,t) — {n — {ni,...,n^) e : Nsi < rn < Nti^i} and Ajv(t) = 
Ajv(0,t). These together with (|2.16p - (|2.18p enable us to represent ^Ar(t) given by 
(jl.2[) in the form 

k I 

(2.19) ^N{t)^Y.^'^N{it)+ ^'-^^W 

1=1 i=k+l 

where for 1 < i < fc, 

(2.20) i,^^{t)^N--''' F,{X{n),X{2n),...,X{tn)) 

TieAjv(t/i) 

and for i > fc + 1, 

(2.21) e,,^(i) = iV-''/2 Y F,{X{q,{n)),...,X{q,{n))). 

neAjv(t) 

Next, we define a matrix {Di_j, 1 < i,j < k) which appears in the limiting covari- 
ances formula in our main result below. For any i,j<k set 

(2.22) A J = f 
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where v is the greatest common divisor of i and j, Cij{u) = ii v does not divide 
ah components of u S Z"^ and for i' ~ i/v, j' ~ j/v, 

(2.23) Cij{vu) = J Fi{xi,...,Xi)Fj{yi,...,yj)l\'"^^j^d^iau{xai',yaj') 

with /io being the diagonal measure, i.e. J f{x,y)d^o{x,y) = J f{x,x)d^{x). 

2.2. Theorem. Suppose that the conditions ^2.4^ ~ [^.ll]) and Assumvtion \2. 1\ hold 
true then each random field (,i^N(t), i = 1,2,...,^ converges in distribution as 
N ~^ OD to a Gaussian random field rji[t) . Moreover, {rii(t),ri2(t), ...,r]e{t)) is an 
£-dimensional Gaussian random field such that rii(t), i <k have covariances 

Er]i{s)Vj{t) = A,j ]^niin(s/,t;), i,j < k 

1=1 

with matrix Di j defined by \2. 2S\) while the random fields rjiit), i > fc + 1 are 
independent of each other and of rjj 's with j < k and have variances 

E\rii{t)\'^ ^ J \Fi{xi,X2, ■■■,Xi)\'^dfi{xi)dfi{x2) ■ ■ ■ dfi{xi), i>k + l. 

Finally, S,N(t) converges in distribution to a Gaussian random field S^{t) which can 
be represented in the form 

k I 

(2.24) ^it) =Y,'n,{it) + ^'W- 

i=l j=fe+l 

In order to understand our assumptions observe that vjq,p is clearly non- 
increasing in q and non-decreasing in p. Hence, for any pair p,q > 1, 

^q,p{n) < il>{n). 

Furthermore, by the real version of the Riesz-Thorin interpolation theorem or the 
Riesz convexity theorem (see [13], Section 9.3 and [llj. Section VI. 10. 11) whenever 
e [0, 1], 1 < po,pi,qo,qi < oo and 

i_i-6i 9 i_i-6i e 
p po Pi q qo qi 

then 

(2.25) m,,p{n) < 2(n7,„,p„(n))i-»(ti7,,^p, (n))^ 

In particular, using the obvious bound zuq-^^p-^^ < 2 valid for any qi > pi we obtain 
from (I2.25P for pairs (oo, 1), (2, 2) and (oo, oo) that for all q > p > 1, 

(2.26) ^^9,p(n) < (2Q;(n))p~9, n7,^p(n) < 2^+p " ? (p(n))^" p + ^ 

and Wq,p{n) < 2^^p{(p{n)y^'p. 
We observe also that by the Holder inequality ioi q > p > 1 and a <E {0,p/q), 

(2.27) /3,(r) < 2i-«[/3p(r)]"7l^ 

p~ga 

with 76/ defined in p.l2p . Thus, we can formulate Assumption 12 . 1 1 in terms of more 
familiar a, p, (j), and i/i-mixing coefficients and with various moment conditions. It 
follows also from (|2.25l) that if zuq^p{n) — as n — > oo for some q> p > 1 then 

(2.28) vjq^p{n) — 7> as n — > oo for all g > p > 1, 
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and so (|2.28|) holds true under Assumption 12. II 

In order to prove Theorem 12.21 we wiU represent £,i^N{t) in the form 
Si<i<Af "where now I is one dimensional which together with estimates of 

the next section will enable us to apply central limit theorems for mixingale arrays 
(see [19] and [20]). This will lead to Gaussian one dimensional distributions in the 
limit but combining this with a kind of the Cramer- Wold argument, covariances 
computation in Section |4] and tightness estimates of Section [5] will yield appro- 
priate Gaussian random fields as asserted in the theorem. Recall (see ^7\), that 
already in the one parameter case z/ = 1 the process S.{t), in general, does not have 
independent increments so also in the random field case ^(t), in general, is not a 
multiparameter Brownian motion. 

2.3. Remark. As a part of tightness estimates of Section [5] we will see that 
sup^>i tg[g 1],. ii^|^i,Ar(t)|'' < C < oo. Hence, applying the Borel-Cantelli lemma 

we obtain as a byproduct that if Si^N = A^'^^^C^w (*) and Sn = A^"^^ Z]i=i Ci,Af (*) 
then with probability one 

lim Si lyit) = for each i, and so lim Siyit) — 0. 

N^oo ' Af-foo A" 

Still, we observe that this strong law of large numbers can be obtained under more 
general circumstances here since, in particular, we do not need for it convergence 
of covariances derived in Section |4] which requires, for instance, more specific as- 
sumptions on (7j's. 

3. Blocks and mixingale type estimates 

We rely on the following result which appears as Corollary 3.6 in [TTI. 

3.1. Proposition. Let Q and % he a-subalgebras on a probability space {fl,J-,P), 
X and Y be d- dimensional random vectors and f = f{x,U!), x G R'' be a collection 
of random variables measurable with respect to % and satisfying 

(3.1) \\f{x,u)~}{y,uj)\\, < C{l + \xY + mx-y\'^ and \\f{x,uj)\\, < C{l + \xn 
where q > 1. Set g{x) = Ef{x,uj). Then 

(3.2) \\E{f{x, ■)\g) - g{x)\u < c(i + \\x\\i+l^^){m,jg ,n) + II A - E{x\g)\\l) 

provided ^ >0, ^ > ^ + + ^ with c = c{C, t, k, 6,p, q, v,d) > depending only 
on parameters in brackets. Moreover, let x = {w, z) and X — {W, Z), where W and 
Z are di and d— di- dimensional random vectors, respectively, and let f{x,uj) = 
f{w,z,uj) satisfy iS. 1\) in x = {w,z). Set g{w) — Ef{w,Z{Lj),Lo). Then 

(3.3) \\E{f{W, Z, ■)\g) - ~g{W)\\. < c(l + \\X\\'^X2)) 

X Kp(e,H) + \\w- E{w\g)\\l + \\z- E{z\n)\\l). 

We will use the following notations 

(3.4) Fi^n,r = Fi^ri.r{xi,X2, ■ ■ .,Xi^i,Uj) = E(Fi{xi,X2, ■ ■ . , , X(n) ) | J"[/^(„) ) , 

Yiiqiin)) = F^{X{qi{n)), . . .,X{qi{n))) and Yi{m) = if m 7^ qi{n) for any n, 

Xr{n) = E{X{n)\J^if^(^n)), Y^^r{qi{n)) = i^i,gi(n),r(-'^r(gi(n)), 

. . . , Xr{qi-i{n)),uj) and Yi^rim) = if m 7^ qiin) for anyn; 
Y{n) = Y,{n) - EY,{n), Y,r{n) = Y^,r{n) - EY,,r{n). 
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For each I e Z+ introduce cubes = {n = {ni,...,ni,) E : {) < rii < 

I for 1 = 1, v} and for / < Z we set also T(/, I) — □(^) \ □(O- Fix some positive 
numbers ^(r + 1) < 877 < t < 1 and set a(l) = 0, 6(1) = 1 and for j > 1, 

(3.5) a(j) = bij - 1) + [(j- - If^l bij) = a(j) + [^1 and r(j) = [f]. 

Set Aj^'(t) = ANit/i) if j < fc and Aj^'(t) = AAr(t) if i > fc + 1. We define now 

(3.6) V^.t,N{l) = E„6A(;'(t)nT(a(i),6(0) ^^.-W^W) 

and W,,t,w(0 = E„6AW(t)nT(6(0,a(i+i)) ^^.KO^H) 

The sets T {b{l) , a{l + 1)) will play the role of gaps between T{a{l),b{l)) and T(a(/ + 
1), 6(Z + 1)) and we will see that the random variables Wi^t,N{l) can be disregarded 
for our purposes while dealing with the random variables Vi^t,Nil) we will take 
advantage of our mixing conditions in order to show that their centered versions 
Vi,t,N{l) = Vi^t,N{l) — EVi^t,N{l) satify mixingale estimates (see [19] and [2^) with 
respect to the nested family of cr-algebras = J-pi^i), I = 0,1, 2, ... where 

r,(/) = {n e Z"; : dist(n,U,<,g,(n(6(0))) < r{l)} 

and we take Q^^^ to be the trivial cr-algebra {0, D,} for I < 0. 
Namely, for any u d N we have 

(3.7) \\E{V.,t,NmlX)h < E„eA„(*)nT(a(0,6(0) mY,^ril){q^{n))\glX)h 

< |T(a(/),6(/))|max„gT(a(0,6(0) P(^',,r(ofe("))l^i-„)ll2 
where \A\ for a set A denotes its cardinality. Next, for u > Z, 

(3.8) i?(y,,,(o(g.(n))|g«J = 

while for all u > 2 by the Cauchy inequality and the contraction property of con- 
ditional expectations 

(3.9) mY,,rMn))\glt)h < 2\\E{Y,,ril){q^{n))\gl%)h. 

Observe that if n G T{a{l),b{l)) then X = (^Xr(^i^{qi{n)), Xr(i){qi-i{n))) is 

Qi -measurable and for large / we obtain also by the definition of qi for i < k 
and by (g^ and (I^TT|) for i > k that 

(3.10) dist(r,(/-2)ur,_i(0,(z,(n)) >{i-iy 

taking into account that a{l) > — 1)^^'^ ^ (^ ^ !)• We can write also that 

(3.11) 

6(0 < (? + 1)1+^, |n(n| = (i + i)''' and |T(f,oi < i^a - i)(i + ir-^ if ? > l 

1 + T 

Thus, applying Proposition 13.11 to the right hand side of p.9p with Q — 
•^ri(i-2)uri_i(i) and H = J^u^^i^(qi{n)) together with (|3.10p . (|3.1ip . Assumption UTT] 
and the contraction property of conditional expectations we obtain that 

(3.12) \\E{Y,^riMn))\gl%)h < ||i?(r.,.(o(g.(n))|e)||2 
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for some C, C > independent of n, I provided p, q and S satisfy the conditions of 
Assumption [211 Collecting (f3Jl) - (|3Jl) . ([3lT|) and ([332]) we obtain that for any 
u>2, 

(3.13) \\Em,t,N mil Jh < C2''(2+-)r(2W2)-i^^_^((/ _ 1)- _ 

In order to incorporate (j3.13l) into the setup of mixingale arrays from |20j we 
consider the triangular array Vi.t^N{l) = A^"''/^V'*,t,Ar(0, I = 1, 2, -L(iV); N = 
1, 2, ... where L{N) = min{j' : a{j + 1) > N}. Observe that 

iV> [f]>il + r)-\LiN)^iy+\ 

l<i<L(JV)-l 

and so 

(3.14) L{N) < {N{1 + t))1/(i+^) + 1 < 27Vi/(i+^) + 1. 

Employing Lemma l4?2] from the next section together with p.lSp . p. lip and p.l4p 
we obtain that 

(3.15) m,t.Nml < C\T{a{l),bm < ^r{l + l)(i+-)('^-i) < CiN'^-^ 

1 + T 

for some C, Ci > independent of I < L{N),N, i and t. Observe that by the choice 
of r and f] we have that z/(2r + ry + 2) — 1 — But < — 1 which together with (|2.13p 
and p. 131) enables us to write 

(3.16) \\E{V,^t,N{l)\Gl'}j\\2 < C^r^ < C^u-^ 

for some C2 > 1 independent of l,N,i,t and u = 2,3,...,L Set -tpj = 
C2(max(l,j))-i for j = 0,1,2,... and m^N = CiiV^^aTTT for I = 1, 2, L(7V). 
Then p.Sp . p.lSp and p. 161) yield the first standard mixingale estimate 

(3.17) \\Ei%t.N{i)\gl'2jh < f^/.wv^u 

for all It = 0, 1, 2, ... and the conditions imposed on ipj^s and tT/,Ar's in [20] can be 
easily verified, as well. The second standard mixingale estimate (see [50]) is trivial 
in our case since Vi^t,N{l) is t/;'^\j-measurable for any u> 1, and so 

(3-18) Em,t^N{i)\gl'lJ - %tM{i) = 0. 

Next, we estimate contribution of small blocks (gaps) Wi^t,N{j), j > 1- Set 

r,{l) = {n e : dist{n,U,<^q,{a{a{l + l)))) <r{l)}, 

G = -^fi_i(0uf i(i-2) and H = ^u,^^,-,(qi{n)) where n g T {b{l) , a{l + I)) . Observe that 
by the properties of g^'s for any such n and large enough I, 

dist(f,_i(0uf,(/-2),(7,(n)) > r. 

Furthermore, if j < Z — 2 then Wi^t,N{j) is J^p.(-;_2')-nieasurable, and so employing 
Proposition [ST] with such Q and H we obtain similarly to p.l2p relying also on the 
Cauchy inequality that 

(3.19) \EW,,t,N{.m^ril){q^{n))\ = |i?(W^,,t,w(j)£;(r,,,(,)(g,(n))|g)) | 
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for some Ci > independent of t, N, n, I and j < I — 2. Hence, by p.6p and p. lip . 

(3.20) \EW,,t^NU)W^J,N{l)\ = \E{W,,t^N{j)E{W,,t,Nil)\g))\ 

< |T(6(0,a(/ + l))|max„eT(fc(0,<^(;+i)) l^(W^M,A^(j)^(>^^KofeW))l^)) I 

< C2||M/.,t,Ar(j)||2Z''(2^+"+2)+2" — ltn,,p(r - l^)) 

for some C2 > independent of t, N, n, I and j < I ~ 2. Now we can write 

(3.21) i?(E2o^ W.,,Mj)f < E'io^ {SEWl^tMj) 

where by (|2.13p . ()3.20p and the choice of r and 77, 

^"^'"+''+'^wp.?(^"-?'') <^?2^/'""tnp,,(r-;'') <oo. 

i<;<L(Af) 1=1 

Relying on (12.31) . (12. 5p . (|2.15l) and the Holder inequality we can estimate the error 
of replacement ofYi{qi{n)) by its r(^)-approximation ^(j-) ((/^(n)) (see Lemma 3.12 
in [17]), 

(3.22) miq^in)) - r,,,(,)(g,(n))||2 < C4/?,^(r(j)) 
for some C4 > independent of i,j and n. Now, set 

1<1<L{N) 

Then by (|mD . (IX^ and (13221), 

(3.23) ||C.,Ar(i) - C..A^(i)l|2 < C^5iV-''/'(Ei<i<L(^) /'^(^+i)-i/3,^(n) 

+(Eo<;<Lw(iiw^M,iv(oiii + m,tMm2)f). 

It follows from ()2.15p and Lemma 14.21 of the next section that 

(3.24) |lVF,,t,w(Z)|l^ = O(|T(6(0,a(? + l))|). 

By AssmTLption l2.1l and the choice of 77 and r we obtain from (|2.14p . (13. lip . p.l4p . 
(15^ and ((X^ that 

(3.25) ||C.,iv(i)-aA'(i)l|2 <C6iV-''/2(Ei<z<L(Ar)^'^^"+'-"'^-' 

< C7(iV''^5- itt) + AT^T+T) ) ^ as ^ CX3, 

and so for each t the limits in distribution as iV — > 00 of Ci.Nit) and of Ci,Af(^) 
coincide (if they exist). 

Observe that similarly to 1^1}, (IXT^ and it follows from f^TT^ . (|XTT|) 

together with the contraction property of conditional expectations that 

(3.26) mMt)\ < N-''^^ Ei<i<LiN) m,t.N{i)\ 

for some Cs,Cg,Cio > independent of N, and so |i?Ci,7v(^)| — > as — ^ cx). 
Hence, the Gaussian limiting behavior of Ci,Jv(0 ~ EQ^N{t) which we will derive 
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via mixingale limit theorems yields the same Gaussian limiting behavior of Ci.Nit), 
and so in view of p.25p the same holds true for ^i.N{t), as well. 

4. Limiting covariances 
The first step in our limiting covariances computations is the following estimate 

of 

bij{m,n) = EYi{qi{m))Yj{qj{n)), TO,n e Z!^ 
where Yi{qi{n)) was defined in (|3.4p . 

4.1. Lemma, (i) For i, j — 1, ...,£ and any m, n G Z"^ set 

(4.1) Sij{m, n) ^iam(inmi<i<j \qi{m) - qi{n)\, \m\) 

and Sij{m,n) — max(sij (TO, n), Sj^i(n,m)). 
Then for all i,j < k, 

(4.2) s^.J{m,n) > ^\im- jn\. 

Furthermore, if i > k + 1 then for any e > there exists Mg > such that if 
max(|m|, |n|) > Af^ and m ^ n then 

(4.3) Si^i{m, n) > min (|m — n\ + e^^ , max(|m|, |rt|)) > — |m — n|. 



(ii) There exists a function h{l) > defined on integers such that X^l^i I h(l) < oo 
and for any i, j — 1, 2, i and 1 — 0,1, 2, 

(4.4) sup \bi,jim,n)\ < h{l). 

Proof, (i) Let i,j<k and set u = \im — jn\. Then i\m\ + j\n\ > u, and so cither 

> ^ or \n\ > jj. Suppose, for instance, that |r7i| > ^. If Sij{m,n) > 
then Si,j{m,n) > and we are done. So assume that Si_j{m,n) < j^. Then 
mini</<j \im — ln\ < j^, and so \im — ln\ < for some I < j. Then i\m\ — l\n\ < 
whence 

1 , . . u u u u 

n > — I m ^ > — ^ > — . 

? ^ ' ' Ak^' - 21 4fc2 - 4k 

Next, let mini</<i \ jn — lm\ = \jn — lm\, I < i. Then 

\jn — lm\ = \ jln — lim + lim — lm\ > [li — l)\m\ — l\ln ~ im\ 

>\m\-k^>f^-^^>^^. 

It follows that Sj^i(n,m) > and so by above either Sij{m,n) > or 
Sj^i{n,m) > whence Sij{n,m) > j^. The case \n\ > ^ is dealt with in 
the same way exchanging i and j as well as m and n in the above argument. 

In order to obtain (|4.3p we rely on the definition (|4.1I) together with the assump- 
tions dm and (HjD). 

(ii) By (|2.10p there exists M such that \qi{m) ~ qi{m)\ > \ni\ for alH < i provided 
|m| > M. Hence, for such m, 

min ( min \q^{m) - qi{n)\, min \qi{m) - qi{m)\) > Sij{m,n). 

l<l<j l<l<i 

Assume that Sij(rn,n) — Si^j{m,n) > 2r and set 

6g(m,7i) = EYAq^imWJ,riqJ{n)) 
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where Fi,r(0 was defined in (fO|) . It follows from (|23| . (|2?5|) . (|2Jl) together with 
the Holder inequality (of. Lemma 3.12 in [17 ) that 

(4.5) |6S(m,n) - h,,{m,n)\ < C{P,{r)f 
where a constant C > does not depend on i,j^m,n and r. Set 

Tr{m,n) = uj,-\C/^(g„(m)) U^„^i i7^(g„(n)) 
where Up{x) = {y : \x — y\ < p}. Applying Proposition 13. II we conclude that 

(4.6) \bi]im,n)\ = |i?(i?(r,,,(g,(m))| J-r„(„,.„))r,,.(gj(n))) | 

< \\Yj^r{qj{n))\\2\\E{Y^,r{qi(m))\J^Y'Ani,n))\\2 < Cr" mq^p{s^^J {m, u) - 2r) . 
The estimate in the case Si.j{m,n) — Sj.i{n,ni) > 2r is similar. Now we choose 
r = |;Sij (m,n) and take 

(4.7) h{l)^C{l-'mp^g{2l) + {l3pJl)f) 

provided either Sij{m,n) = Sij{m,n) and \m\ > M or Si,j(m,n) — Sj,i(n,m) 
and \n\ > M. Next, suppose that Sij{m,n) = Sij{m,n) and |m| < M. Then 
Sj^i{n,ra) < Sij{m,n) < M, and so either \n\ < Mormini<i<i \qj{n) — qi{m)\ < M. 
In the latter case we must have \n\ < N{M) for some N{M) > depending only 
on M. The same argument holds true when Si j{m,n) — Sj^i{n,m) and |n| < M. 
Thus, when I > M we can define h{l) by (|4.7p while for Z < M we can take 

h{l) = max{|6ij(m,n)| : < |m|,|n| < max(M, iV(M)), I < ij < i} 

concluding the proof of the lemma. □ 

For s,te [0, 1]" with s < t set 

C,,^(s,t)=7V-'^/2 ^ F,(X(<zi(n)),...,X(g,(n))). 

Now we can obtain an appropriate estimate of the second moment of ^i,N- 

4.2. Lemma. There exists C > such that for all t = (ti, t,^) > s — (si, s,^) > 

and i — 1^ ^, 

(4.8) i?|6,iv(s,t)P < CiV"nAA,(s,t)| < C J](ii - Si +7V-1). 

(=1 

Proo/. Set G(m, = {r^ e Z+ : / < |m - n| < ^ + 1}. Then by (g^l) and gH) for 

1 < fc, 

(4.9) i?|6,jv(s, OP = ^-"^ E™,„eA„(.,t) &m(™, ") 

< 2iV~-Ez=oM[4S-])E™eA„(M) |G("^'OI- 
If i > A; then by (g^l) and 

(4.10) ii;|c,,^(s, i)p < TV-"- E„,„:|™i>i<M, l&^»("^. ")l 

+2^"" E^O Em6A„(s,t) E„GG(m,;) I^V^C"^' I 

< N--{C^ + E;:o ^^([^2]) E™eA„(.,t) |G(™, 01) 
for some Ci > independent of s, i and N . Clearly, for any / > 0, 

(4.11) \G{m,l)\<C2(l^\r-^ 



12 



Yu.Kifcr 



for some C2 > independent of m and I since it is bounded by the number of 
integer points between spheres of radius I and I + 1. Finahy, by (|4.9|) - (|4.11|) and 
the summability of l'^~^h{l) we derive (|4.8p . □ 



4.3. Lemma. For each t > and i,j < k, 

L 



hm HmsupTV"'' V \b, An,n')\ ^ Q. 

U^OD JV_).oo 



Proof. Set G,,j (n, I) = {n' E Z"^ : I < \in ~ jn'\ < / + 1}. By LcmmaSH 
(4.12) Ar-''Eo<...'<«t 

<n<Nt 

for some constant C3 > independent of u and A'^. As u — > 00 the right hand side 
of (|4.12p tends to zero in view of summabihty oi l'^ h{l) in I . □ 

4.4. Proposition. For any i,j < k and s — {si, s„),t = {ti, ...,t^) > the limit 

s-^ , ,^ vXY, ,min(s/,t;) , , 

Q<in<Ns *^ 
0<jn'<Nt, in-jn'=u 

exists for any u € Z"^ where v is the greatest common divisor of i and j , Ci,j{u) = 
if V does not divide all components of u E and Cij{vu) was defined by h2.2S\) . 
Finally, 

(4.14) hmjv->oo £'6,Af(s)0-Jv(*) = limAf^oo N^" J2 o<."<«= bij{n,n') 

Q<jn' <Nt 

= Ajnr=i™"^(*''*') 

with Di^j defined by i2.22\) . 

Proof. Let v be the greatest common divisor of i and j. If w e has components 
which are not divisible by v then the sum in (j4.13p is empty, and so in this case 
Cij{u) = 0. Thus it remains to deal with this sum when in — jn' ~ vu for u € . 
We will show first that the hmit 

(4.15) Ci^j(u) — hm bij{n,n') 

\n\,\n'\^QO.^n—jn'—^}u 

exists. Observe that if we consider two strings {n,2n . . . , in) and {n', 2n', . . . ,jn') 
with in — jn' = vu then there will also be pairs (zaj Ja), a = 1, 2, . . . , f — 1 such 
that ian — ja n' — au where ia — aii and ja — aji with ii and ji being coprime. 
On the other hand, if i/j ^ i\l j\ then 

(4.16) |in-jn'| ^ j| in - n'| - j| (j - ^)n + ^(iin - jm')] 

> Irtl?!- — -y-l — ^lul — > 00 as \n\ — > 00. 

We split the collection of numbers (n, 2n, m; n', 2n', ...,jn') into disjoint sets 
Fi, Ftj, Fi+j-i, where Fq, = {aii,aji}, a = l,...,v are pairs and /3 — 

l,2,...,i + j — V are singeltons. We order the latter so that F^ = {Z^n}, 1 < 
1/3 < i, P — V + l,...,i with Ip 7^ ail for a — l,...,v and T^+p = {I' pn'} for 
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P = i + 1, i + j — V, 1 < I' fs < i, I'p ^ aji for a ~ 1, v. By ()4.16p there is 
(5 > depending on u but not on n and n' such that 

min distfF;, F;/) > 

Set 

[/^(r,) = {n e l.\ : dist(n,r;) < r}, / = l,2,...,i + j - u 

and choose r — r(n) — > oo as — )■ cx) so that all Ur(n)^l)^ I — Ij 2, i + j — u 
were disjoint. 

Now, observe that bi,j{n,n') has the form 

EG(Vi{n,n'),Y2{n,n'), ...,Yi+j^^{n,n'j) where Ya{n,n') = {X{ai,n),X{aj,n')) 
for a = 1, V, Yfi{n, n') — X{lpn) ioi P ^ v + 1, i and Yp{n, n') — X{l'pn') for 
/3 = i + 1, i+ j ~ V. Define Gi = G and successively 

Gi+i{yi +1 1 7 • ■ • 7 yi-\-j—v ) = EG{Yi{n, n'), yi+i,yi+2, Vi+j—v)- 

Relying on the assumptions (|2.13p and (|2.14p we can apply p.3p of Proposition [SH] 
with V = {Yi+i{n,n'),Yi+2{n,n'), ...,Yi+j^^{n,n')), Z = Yi{n,n'), Q ^ with 
U = ^i+i<p<i+j^vUr{n){^ p) and % = -^(7,(r(n)) obtaining that 

EGi{Yi{n,n'),Yi+i{n,n'),...,Y,+j^^{n,n')) 
~EGi+i{Yi+i{Yi+i{n, n'), n')) as rn- cx). 

This argument repeated for I ^ 1,2, i + j — v — 1 yields ()4.15p with c^.j given by 

Finally, in order to obtain (j4.13p and ()4.14p we have to count the number of 
solutions n, n' of the vector Diophantine equation in — jn' — vu where < in < Ns 
and < jn' < Nt. Since we have to satisfy this equation coordinate wise, the 
number of solutions is the product of the number of solutions in each coordinate. 
Let, as before, i — vii and j = vji with ii and ji being coprime then all solutions 
of the equation iini — jin'i = ui are given by n = no+jim and n' — nQ + iim where 
no, n'fj is its particular solution and m is any integer. The number of such solutions 
with < in; < Nsi and < jn'i < Nti for large iV is equal approximately to 

AT w *^^_Lnn^ Nvmm{si,ti) 

Jy min( — , -T^) + 0(l) — h (J[l) 

and taking the product in I we obtain (I4.13P while (|4.14p follows from (|4.13l) and 
Lemmas WJ\ and □ 

4.5. Proposition. For i > k + 1, 

(4.17) lim = (TTt;) / (F.i{xi,X2, ■ ■ ■ ,Xij) dn{xi)dfi{x2) ■ ■ ■ dfi{xi). 

1 = 1 '' 

Moreover, for any t, s £ and j < i, 

(4.18) lim ^(e,,Ar(t)0,iv(s)) =0. 

Proof. By (l^- dmi) . and gll, 

bi^i{m, n) — > as max(|m|, |n|) — > oo with |m — n| > 1. 
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Therefore, by ()4.3p and ()4.4p for any fixed L, 

liniSUp^_,^iV-'^X;m,„6A„(t),rn#nl^^^("^''^)l 

< limsup^^^iV-''(X;i<|™_„|<L \Ki{'m,n)\ 

+c\^N{t)\Y.i>Li''-'m)^cY.i>L/2^''-'Ki) 

for some C, C > independent of N and L. We now iet L — >■ oo and since l'^~^h{l) 
is summable it follows that lim sup in the left hand side above equals zero, i.e. the 
off-diagonal terms do not contribute in (|4.17p . It remains to deal with the diagonal 
terms 6i,i(n, n). Since \qj{n) — gj_i(n)| oo for j = 2, 3, ...,£ as \n\ — >■ oo it follows 
by the argument similar to one applied in the proof of Proposition 14.41 (see a more 
general Lemma 4.3 in [IT]) that 

(4.19) lim 6j_i(n,n)= {Ft{xi,X2, ■■.,Xi))^d^{xi) ■ ■ ■ dn{xi). 

n|-i.oo J 

Namely, set Gi{xi, ...,Xi) = (^Fi{xi,X2, ■■■,Xi))'^ and recursively for Z = i — 
1, ...,2,1,0, 

Gi{xi,...,xi) = J Gi+i{xi,...,xi+i)dfi{xi+i). 

Taking into account that \qi{n) — qi{n)\ — > oo as |n| — >■ cx) when I ^ I we apply 
Proposition 13 . 1 1 to obtain successively for Z = i, i — 1, 1, that 

|£;G,+i(X(qi(n)),...,X(g,(«)))-SGz+i(X(gi(n)),...,X(q,(n)))| -^Oas |n| ^ oo. 

Since bi^i{n,n) = Gq we arrive at (|4.19p . 

Next, we deal with (|4TT8l) . Since i > j and i > k then by ([2Jl) - ([2TT|) for any 
£ > there exists N{e) such that whenever |m| > eN and \n\ < Ny/v we have 
Sij{m,n) > min(|m — n| + s~^, |m|) provided N > N{s) where Sij('m,n) is the 
same as in Lemma 14.11 It follows from Lemmas 14.11 and 14.21 that 

(4.20) mMmj,N{s)\ < <6N, ne Ajv (s) ^iJ {m,n)\ 

+N ''Z]|m|>£A',meA„(t),neA„(s) 

for some C > independent of N and e. Letting in (|4.20p . first TV — >■ oo and then 
£ — we arrive at (|4.18p . □ 



5. Tightness estimates 

First, we will extend the estimate of Lemma [4.21 to the corresponding estimate 
of the 4th moment. 

5.1. Lemma. There exists C > such that for all t = (ii, tv) > s = (si, s^) > 
and i — 1, 

(5.1) E%,^{s,t)t < C7V-2-|Ajv(s,t)P < Cl[{ti - si+N-')\ 

1=1 
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Proof. For n^^), n^^), n^^), G Z"; set 
and for r > 0, 

n'^'^') = i?>^i,g,(„(i)),r5^i,gi(n(2)),r^j,g,(n(3)),r^j,9i(nW),r- 

Then similarly to ()4.5p . 

(5.2) |d|'-)(n«,n(2),n(3),n(4))-d,(n(i),n(2),n(3),n(4))| <Ci(/3,(r))* 

where Ci > does not depend on ri'^'^\ n^"') and r. Define 

Wi(r;.'^"'"', n*^^', n*^^', n*^*') = maxi<j<4 ( min ( min;<i \qiin^-'^) — 

minj^^.,<J(?.(n(^))-gz(n(^^))|)). 

Without loss of generality assume that 
(5.3) 

v^(n(i\n(2),n(3),„(4)) ^ min (min - g,(n(^'))|, min |q,(n(^)) - g,(n(j'))|). 

For each a > introduce the sets 

= e Z+ : a < ^;,(n(i), n^^), n^^), n^) < a + 1} 

and 

r47V,s,t) = {(nW,n(2),n(3),„(4)) ^ r,, : n^'\ n^'\ n<-'\ n'-^^ e AAr(s,i)}. 
If i < fc and (n(i),n(2),n(3),n(4)) g then for j = 2,3,4, 

either jn^-'-') < a + 1 or jm'-'-' — Zn'-'-'l < a + 1 for some I = 1, and j 7^ 
It follows that 

(5.4) \raiN,s,t)\<C2a^{l+a^ + \AN{s,t)\'') 

for some C2 > independent of a, N, s and t. If i > fc + 1 then by ([2Jl) - ([2TT|) 
there exists M > such that whenever \n\ > M, 

(5.5) min \qi{n) — qi{n)\ > \n\ and min \qi{n) ~ qi{h)\ > \n — n\. 

Ki i<i, 

Then similarly to the case i < k we conclude from (|5.5p that (|5.4p holds true also 
when i > k + 1. 

Next, let r^a/3 and (n^^), n^^), n^^), n^) e Tq satisfy ([531). Set 

vI/,(n(i),n(2),n(3),nW) = ur=\C/.(g«("^'^)) U ( U.t^ UUt^.(9.("*^'^))- 
Then by Proposition 13.11 similarly to (14. 6p we derive that 

(5.6) |d|'-)(n(i),n(2),n(3),„(4))| < |ii;(ii;(K,^^(„a,)^,| J-^^(„a),„(.),„(3),„(4))) 

><^j,g,(n(2)),r^i,g,(n(3)),r^i,?i(nW),r) I < C's^'^^^^.p (a/3) 

for some C3 > independent of a, n^^^ n'^^ and n^4)_ ggt g((j) — 
Ciipgia/3)y + C3a''zUg,p{a/3). Then 

(5.7) i?|^,^(s,t)|4 < 2- „,3,,„(.,),rj^,M) 

<C2Er=o9(a)a'(l + «' + |A^(s,t)|2) 

and (|5.ip follows from (|5.7p and Assumption 12.11 □ 
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Now tightness of each sequence of random fields {£,i,N{t), t e [0, 1]''} follows by 
a slight modification of |9j (see also Ch.5 in [8] and Theorem 1.4.7 in [16]), and so 
the sequence of random fields {^N{t), t e [0, 1]"} is tight, as well. 

6. Gaussian limits 

For each fixed t e [0, 1]" the convergence in distribution as iV — > oo of each 
Ci,Nit), i = !,...,£ to corresponding Gaussian (maybe degenerated) random vari- 
ables follows from central limit theorems for mixingale arrays (see [5D], [Tl] and 
references there) in view of the mixingale estimates of Section [Sj convergence of co- 
variances obtained in Section |4] and the tightness result derived in Section [5j Then 
£.i,N{t)j * = I7 •■•j^ also converge in distribution to the same Gaussian random vari- 
ables in view of (1X251) . Furthermore, for any t = {t^^\ ■•■,i'-'^), t^"^ e [0, l]"", a = 
1 , . . . , J and d — {di, ...,dj) set 

j 

(6.1) l^^,t.d,iv(0=E^-^M(").w(0 

a=l 

and 

(6.2) W(t) = ^d,C.,Jv(i'"))-iV-'^/2 ^ V,,t,d,N{l). 

0=1 l<'<i(A') 

Then, we obtain from (j3.13p and (j3.18l) similar mixingale estimates also for 
Vi,t,d,Ar(0 which via [20] and [H] yields convergence in distribution as — 00 
to Gaussian random variables of each Ci,d,jv(t). This together with p.25p imply 
that each 

j 

(6.3) e.d,iv(t) = 5]d,6,Ar(t('^)) 

a=l 

converges in distribution as — >■ cx) to Gaussian random variables. Hence, finite 
dimensional distributions of each ^j^^v have Gaussian limits which together with 
tightness results of Section [5] yields that each ^i^^r converges in distribution as 
TV — >■ 00 to a Gaussian random field rji. 

In fact, we can show that (fi,jvi ^i,k) converges in distribution as -> 00 to a 
fc-dimensional Gaussian random field (rji, 77^). Indeed, for any e — (ei, Cfe) € 
R*^ set 

fc fc 

(6.4) Vt,d,e,Ar(0 = XI ei^j,t,d,A'(0 and Cd,e,jv(t) = X eiCvd,A'(t). 

i=l i=l 

Then it is easy to see again by (j3.13p and p.lSp that similar mixingale estimates 
hold true also for Vt,d,e,jv(0 which via [20] and [14] provides convergence in distri- 
bution as iV — >■ 00 of Cd.e.iv to a Gaussian random variable which must have the 
same distribution as X^iLi E1=i e.i-qi{t^°-^) . As above we conclude from p.25p and 
tightness arguments of Section[S]that, in fact, X^iLi ^i^i,N converges in distribution 
as ^ cx) to a Gaussian random field which must have the same distribution as 
Si^i- Thus, (771, T^fe) is a Gaussian random field and (^i.at, Cfe,Af) con- 
verges in distribution to it as A^ — > 00. Finally, X^iLi converges in distri- 
bution as A^ 00 to the random field X]i=i Vii'^t) which must be Gaussian as a 
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result of the linear transformation (in the path space) of a Gaussian random field 
(see, for instance, |5, Section 2.2). 

Next, clearly, ^a? converges in distribution to ^ given by (j2.24p and it remains 
to show that rji with i > k + I are independent of each other and of rji with i < k 
which will imply that ^ is a Gaussian random field. This can be done either via a 
modified version of Theorem 5.6 from [17, or by the following more direct approach. 
First, observe that ()2.11|) implies that there exists ejv — as iV — cx) such that 

(6.5) lim mill (|g,(n)| - max |9,(n)| - |n|) = oo 

Af-s-oc n,neAN(l)\AN(ejvl) 

where 1 — (1, 1) g [0, 1]". We see from (|6.5p that if n,n e A(£„l,t) then any 
qi{n), i > fc + 1 and qj{n), j ^ i are widely separated. Next, it follows from Lemma 
14.21 that for each i, 



(6.6) ^i,Af(£Afl) — in probability as — > oo, 

and so in all our arguments the sum over Ajv(£Arl) can be disregarded. 

In order to use (|6.5p and (|6.6p introduce jn — min{j : a(j) > e^N} with a{j) 
defined in (03]). For t = (i^^), t^^)), t^-^) e [enAY, a = d = {di,...,dj) 

and e = {ck+i, ei) define for 1 — 1,2, L{N), 

3 I 

(6.7) Vt,d,e,iv(0 = E E dae.V,^,(a,^^{l) 

a=l i=k+l 

as in dm and (H^]). On the other hand, for {i-k)L{N)+jN < I < {i + l-k)L{N) 
and i > + 1 we set 

i 

(6.8) Vt,d,e,jv(0 = ^ dae,F,_4(a,^^(Z - (^ - fc)L(iV) 

setting 14,d.e,Ar(0 to be zero if it is not defined by one of the above and assuming 
that > e^v for aU a. It is easy to see from ([^ - (PH]) . (15^ and mixingale 

estimates of Section [3] that the sequence (Vtde7v(0 ~ EVt a e n{1))N~^/'^ , I — 



Jn, ■■■,iL{N) forms a mixingale array with the corresponding cr-algebras Qi — 

for i < k, Gi+{i-k)L{N) — Gi^^ for i > k + 1 and estimates similar to p.l7p . Obsen 
that 



Cd.e,Ar(t)-iV'^/2 ^ Y^e,V,^t,d,N{l)^N-'/^ ^ Vi,d,e,iv(0- 

jN<l<L{N)i=l 1<1<IL(N) 

Taking into account (|6.6p , mixingale estimates of Section [3] and convergence of 
covariances results of Section |4] and tightness arguments of Section [5] we conclude 
by mixingale limit theorems from [2D] and [13] that Cd.e,Af(t) — -ECd,e,Af(t) converges 
in distribution as A'^ — >■ oo to a Gaussian random variable while similarly to the end 
of Section|3]we see that Cd,e,Af(t) converges to the same limit. It follows from here 
in view of (I3.25P and (|6.6I) that each linear combination X^^^i Y^\=i daSi^i^Nit'-"'^) 
converges in distribution as — )■ oo to some Gaussian random variable which then 

must be Ea=iELi^^ae.?7.(i'"^)- Thus {m ,7^2, •■•,%) is an ^-dimensional Gaussian 
random field. Invoking again the linear transformation argument from Section 
2.2 of [3] we conclude that the field £^{t) given by (|2.24p is a Gaussian one and 
taking into account the vanishing limiting covariances result (|4.18p we obtain also 
independency of ?7i, i > fc + 1 of each other and of rji^s with i < k. □ 
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